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Surface mass density of the Einasto family of dark matter haloes: 
Are they Sersic-like? 
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ABSTRACT 

Recent advances in N-body simulations of dark matter haloes have shown that three-parameter 
models, in particular the Einasto profile characterized by dlnp(r)/dlnr oc with a shape 
parameter a < 0.3, are able to produce better fits to the 3D spatial density profiles than 
two-parameter models like the Navarro, Frenk and White (NFW), and Moore et al. profiles. 

In this paper, we present for the first time an analytically motivated form for the 2D 
surface mass density of the Einasto family of dark matter haloes, in terms of the 3D spatial 
density parameters for a wide range of the shape parameter 0.1 ^ a ^ 1. Our model describes 
a projected (2D) Einasto profile remarkably well between and (3 — 5) r2oo, with errors less 
than 0.3 per cent for a ^ 0.3 and less than 2 per cent for a as large as 1 . This model (in 2D) can 
thus be used to fit strong and weak lensing observations of galaxies and clusters whose total 
spatial(3D) density distributions are believed to be Einasto-like. Further, given the dependence 
of our model on the 3D parameters, one can reliably estimate structural parameters of the 
spatial (3D) density from 2D observations. 

We also consider a Sersic-like parametrization for the above family of projected Einasto 
profiles and observe that fits with a Sersic profile are sensitive to whether one fits the pro- 
jected density in linear scale or logarithmic scale and yield widely varying results. Structural 
parameters of Einasto-like systems, inferred from fits with a Sersic profile, should be used 
with caution. 

Key words: gravitational lensing - galaxies: clusters: general - galaxies: fundamental pa- 
rameters - galaxies: haloes - galaxies: structure - dark matter. 



1 INTRODUCTION 

Gravitational lensing signatures are a response to the projected sur- 
face mass column density of matter E(_R) along the line of sight 
in galaxies and clusters. Upon a suitable deprojection and circular 
averaging, an estimate of S(7?) can, in principle, be used to trace 
the spherically averaged 3D density profile p{r). 

In the past few years, N-body simulations have shown 
Power et alj f2003). Navarro et al. (2004) (Nav04). iMerritt et al.l 
2006h (M06). IStadel et.al.1 (i2009)(S09)] that three-parameter 



models, especially t he Ei nasto CEinasto. .1965 ) profile and the 
IPrugniel and Simiei] ( ll997h de-proiected Sersic profiles (PS97), are 
able to produce better fits to the 3D density profiles of galaxy 
and cluster-sized dark mat t er hal oes than two-para meter models 
I Navarro. Frenk and Whit3 ( Il997h & (1996) (NEW). iMoore et al.l 
1 1999h ). While th e PS97 profil e has a well known 2D sky projected 
form - the Sersic jSersid 19681) profile, there has been no such ana- 
lytical counterpart for the Einasto profile. 

In this paper, we present a very good approximation for the 2D 
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projection of the Einasto family of 3D profiles. Thus, if the 3D total 
mass density is believed to be Einasto-like, our model can be used 
to parametrically describe the projected 2D surface mass densities 
of galaxies and clusters in the weak and strong lensing regimes. 
However, note that even upon radial averaging to smooth out the 
substructure, not all haloes subscribe to an Einasto profile. Eor the 
rest of this paper, we will limit the discussion to the 2D projection 

of the special case where the 3D profile is Eina sto-like. 

In 3D, the functional form of the Einasto teinastd[T96^ pro- 
file is given by: 



ps Ts 



(1.1) 



where, p{r) is the 3D (spatial) density at r, n (or a = —) is the 
shape parameter, 6 is a function of n, ps the spatial density at a 
scale radius and p(0) — psc'^^ 

In 2D, the Sersic (Sersic 1 9681) profile, which has been used 
to describe the projected surface brightness profiles of galaxies, is 
similar in form to the Einasto profile in 3D l ll.lt and is given by: 



Inf 



1] 



(1.2) 



such that Eg(0)=E_R£, where, R is the projected distance in 
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the plane of the sky, Eu^ is the Une of sight projected surface 
brightness at a projected scale radius Re, which can be defined 
to be the half-light radius of a Sersic profile under the condition: 
q = 2m- 0.3333 + 0.009876/m (PS97) with m (or A=i) char- 
acterizing the shape of the Sersic profile. 

The parameter a—^ in Jl.lb defines the shape of the Einasto 
profile. In the first ever fits to N-body haloes with a Einasto profile, 
Nav04 found an average value of 0.172±0.032 for a wide range 
of halo ma s ses fro m dwarfs to clusters. For galaxy-sized haloes 
IPrada etalJ ( l2006h found 0.133s;a^0.167. iHavashi and White! 
( l2008 h observed an evolution of a with mass and redshi ft (2) in the 
Millennium Simulation {M S) of ISpringel et.alj ( l2005h and found 
Q~0.17 for galaxy and ~ 0.23 for cluste r-sized haloes. A simi- 
lar trend is supported by iGao et al.l ( 12008') where a~0.3 for the 
most massive clusters in MS. Hence, although in this paper we dis- 
cuss our approximation to a projected Einasto profile for O.l^Ja^l, 
particular attention is drawn to the domain a<0.3; where as we 
shall show in §3, the errors due to our approximation out to 30r_2 
(~(3 - 5)r2oo) are < 0.3%. 

The parameters r^, b and n in the Einasto profile are not in- 
dependent. It can be seen, that in terms of a dimensionless length 
X = jr, the logarithmic slope of the density profile is given by: 



13- 



dln{p/ps 
dlnX 



-X 



l/r. 



(1.3) 



Nav04 chooses to define such that b/n = 2 (the isothermal value 
of /3 = —2) at r = and hence label as r_2 and ps as p_2. 
Another approach is to use the convention of M06, requiring 
to include half the total mass. They quote a numerical estimate of 
fe=3n - 0.333 + 0.0079/n. In this paper we will follow the Nav04 
parametrization. 

Since the Sersic profile describes the 2D surface brightness 
profiles of galaxies reasonably well, a natural question is: does 
the 2D surface mass density of Einasto-like 3D dark matter haloes 
also follow a Sersi c-like description? For the Nav04 simulations, 
iMerritt et"^ ( l2005h (M05) have shown that a Sersic function does 
produce fits with acceptable errors in the range of rconv to r2oo, 
where rconv is the minimum radius of convergence and r2oo is the 
virial radius in the N-body simulations of P03 and Nav04. 

In this paper, we focus on the analytical description of the pro- 
jected Einasto profile. In order to see how well a projected Einasto 
profile is described by a Sersic profile, we numerically project i ll. Il l 
and find that a Sersic profile produces acceptable fits only in a lim- 
ited range of the projected radius R. We shall show in §3 that this 
range also depends on whether one fits the numerically projected 
Einasto profile in linear scale (large errors with increasing R) or 
log scale (large errors with decreasing R) of density. It is thus evi- 
dent that the Sersic profile does not give an adequate description of 
a projected Einasto profile for all R. 

The choice of Sersic fit (log or linear) may, for example, have 
possibly strong implications in the strong and weak lensing regimes 
respectively, yielding incorrect results. It is with this perspective, 
and the observation that the projection of an Einasto profile (i.e. 
integral of a Sersic-like function) is not a Sersic-like function but 
rather a Gamma-like function, we present a derivation of the sur- 
face mass density (§2.1) for the Einasto profile. 

Throughout this paper we will describe 3D spatial parameters 
as r, Ts (or r_2), Ps (or p-2), b and n (or a = l/n))and 2D pro- 
jected parameters as R, Re, q and m (or A = 1/m). We will also 
refer to the Sersic profile as Eg, our approximation to a projected 
2D Einasto profile as E_b and a numerically projected Einasto pro- 
file as Eat. Further, for this paper, rconv and r2oo have no physical 



meaning as such. We note from the Nav04 simulations, that on an 
average rconv ~ 5 per cent of r_2, and 7-200 ~ 8 — 10 per cent of 
r_2 for galaxies and ~ 4 per cent of r_2 for clusters. In this paper 
we will, for the sake of discussion, refer to rconv ~ 0.05 r_2 and 
r2oo as ~ 6 r_2. 

In §2 we derive a semi-analytical form ( |2.10| l for E_b, which 
for the Nav04 parametrization of 6 = 2n, simplifies to ( 12. lit . In §3 
we discuss details of estimating the fit parameters in our expression 
for Ee, followed by a discussion of errors in our approximation. 
We then present a comparison between Ejv and the best-fitting Ser- 
sic profile Eg. In §3.4, we note the remarkable accuracy with which 
( 12. lit , along with l |2.12| l and l |2. 13l l, can be used to extract the 3D 
parameters of Einasto-like systems from its 2D form. 



2 SURFACE MASS DENSITY OF THE EINASTO 
PROFILE 

The Einasto profile has generated considerable interest of late and 
has been used in recent studies of dark matter haloes (Nav04, 
Pr06, M06 , G08, HW08, S09). For the N-body haloes in Nav04, 
iMerritt eTal.. (2005) have shown that within the limited radial range 
of rconv ^ r ^ r2oo and shape parameter 0.12 ^ a ^ 0.22, 
a deprojected Sersic profile also fits the 3D distributions almost 
as well as the Einasto profile, and a Sersic profile fits the non- 
parametrically estimated 2D surface densities of dark matter haloes 
with acceptable errors (~ 5 per cent). 

In the following discussion, we argue that if a 3D distribution 
is Einasto-like, the 2D distribution need not be Sersic-like and pro- 
vide an analytically motivated functional form for the 2D projec- 
tion, which can be used to describe the surface mass density of the 
Einasto family of dark matter haloes subscribing to a wide range of 
the shape parameter (1 ^5 n ^ 10 or 0.1 ^ q ^ 1) and over a wide 
radial (projected) range ^ R ^ 6 r2oo. 



2.1 General shape of 2D Einasto profile: A short discussion 

In terms of the line of sight distance (2), the surface mass density 
E(R) can be estimated from: 



E{R) = 2 / p{V + R'')dz 



(2.1) 



While an exact analytical expression for the integral in ( 12. It for 
the Einasto profile dl.ll l has so far eluded us, we derive below an 
excellent semi-analytical approximation. 

To intuitively motivate the functional form of T,e{R) for the 
Einasto profile dLU . observe that at 7? = 0, the integral l |2.1| l 
presents us with an exact solution: 



^£^(0) = ^ r[n] 

where, r[n] is the complete Gamma Function. 



(2.2) 



Hence, for R>0, it is reasonable to expect the integral to depend on 
terms involving incomplete gamma functions. In fact, for the sake 
of discussion, one can make a very crude assumption that most of 
the contribution to the integral in ( I2.lt at a given R (especially for 
R < r-2) comes from the region z > R. Integrating, from some 
i^R to 00 (where > 1), one gets: 



^(^) = ^ r[n,6( — )"] 



(2.3) 
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Similarly, the integral for R > r_2 will have dominant contribu- 
tions from terms involving 7[a,a::]; where r[a,a::] and 7[a,2;] are 
the upper-incomplete and lower-incomplete gamma functions re- 
spectively. 

This is quite unlike a Sersic l ll.2t function. Hence, although a 
Sersic profile may fit a projected Einasto profile in a limited range 
of R it need not be a very good fit for all R. 



2.2 



Sb(R): An analytical approximation of projected 
Einasto profile 



With the 3D spatial distance r = \'z^ + R?, where R is the 2D 
projected distance in the plane of the sky, and z the line of sight 
distance from the object to the observer, at any given R, one can 
define 3 regions for the integral in J2.1b for the Einasto profile 
(ED: 

Region I: z < R, integrating from z = to ^i-R, with ("i ^ 1 
Region II: z > R, integrating from z=^2R to oo, with (^2 > 1, and 
Region III: z ~ _R, in a neighborhood 5 between C^iR and (^2R 

In Region I: z < R, the first term on the right hand side 
(RHS) of JI.ll l can be written as: 



z^+R^ 



R 



r) 



(2.4) 



Neglecting 4th and higher order terms in (z/R), in the binomial 
expansion of i2.4l . the integral of i2A\ (from 2; = to CiR) has an 
analytical approximation: 



psrs exp 



-b 



-1 X =^ 



2n f R\^^'^'' 
b'WJ 




(2.5) 



In Region II: z > R, a binomial expansion of the first term of the 
RHS of (EB gives us: 



2 ,i 



-K-)-(i + (-r)- = ) 



-bi 



— (-)'" 



1 / I .,--/?,4_ 



+ . 



(2.6) 



Fits to N-body simulations with the Einasto profile have so far indi- 
cated an n>3.0 (q<0.3). Hence, observing that the leading contri- 
bution comes from the 1st term b(z/rs)^ , one can drop the remain- 
ing terms. This is especially true in our primary domain of interest 
(3^n^ 10 , i.e. a^CO.3). Even for 1 3, although the approx- 

imation is not as good (as it is for rz^3), it continues to provide 
better fits than a Sersic profile in the entire range l^n^lO. With 
this understanding, the integral from z — C,2R to co, in Region II 
can be written as: 



2e TsUps 



n, b 



C2R 



(2.7) 



In Region III: (^iR ^ z ^ C,2R, one can not make the approx- 
imations made in regions I and II. However, there exists a point 
eR between C^\R and C,2R, where the mean- value approximation 
will be valid in a domain 5 about eR. Since the density profile falls 
rapidly for R « r_2 and gradually for R >> r_2, the domain 
of applicability of the mean- value approximation will be such that 
S < (C2 - Ci ) for 7? < r_2 (-5> as 7? -s- 0) and tending to ^2 - Ci 
for R >> r_2. Further, it should be obvious that ("1, (2 and S will 



depend on the shape parameter n as well. A function describing S, 
with such a property is: 



5 = (C2-Ci)[l-exp(-(ii/r,)'')] 



(2.8) 



with fi = ^i(n) 



The remaining (small) excluded region does not add significantly 
to the integral (refer to error plots in §3). 

With this, the contribution to the integral in J2.I| (, from region 
III- a domain 5 around R - from an application of the mean- value 
theorem at eR can be written as: 



25Rp{z = eR) = 2SRps exp -6 



(2.9) 



A few important observations are in order. First, neglecting 
terms in the integrand of region I l l2.4t and II l l2.6t , leads to over- 
estimating the integrands in those regions. Second, ignoring the 
contribution from region III and fitting only for and 1^2, produces 
good fits (refer to discussion following i2.13t ) but understandably 
with a ^2 greater than (^2 with region III included, i.e. including 
region III, lowers 1^2 allowing it to be closer to 1, resulting in a 
larger contribution from the upper-incomplete gamma function. A 
simultaneous fit of (^1,(2, e and p(m S) accounts for these excess 
contributions through a negative sign from region III. 

In region I, since 4*'' and higher order terms in z/R are ne- 
glected, one can fix ^1 = 1. Further, although e should in principle 
be estimated, we found it to be a reasonable approximation to fix 
e = ^^-j^. This reduces the number of parameters to fit to only 
two - (^2 and /i - which in turn, during the fitting process, compen- 
sates for the approximations made on i^i and e. 

With these approximations, the surface mass density E_e(R), 
for 1 ^ n ^5 10 with X — R/vs can be written as: 



^e{R) 



r(n + l) 



nr [n,6(C2^)*] 




r^'2^^ 



(2.10) 



where. 



Eb(0) 



2e^rsPsV{n + l) 



For the Nav04 parametrization of ^ = 1 and by choosing 
^1 = 1, 1 I2.IOI 1 simplifies to l |2.Il| (: where, after substituting factors 
like6/2n with 1 in l lOOt . we leave b in the rest of the equation in 
order to reduce clutter. This gives us: 



Sb(7?) = 



r(n + i) 



n F 



2' 




with. 



5 = {<:2 -<:.)[!- e-Wrsr] 
, C2 + Ci 
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2 4 6 8 10 12 

n 

Figure 1. C2(n) ED for 0.1 < a 1.0 (1.0 ^ n 10.0). Each point 
has been obtained by fitting the numerically projected Einasto profile with 
tzTll for h = 2n and Ci = 1. 




2 4 6 8 10 12 



n 

Figure 2. n{n) tTTst 0.1 ^ a ^ 1.0 (1.0 ^ n ^ 10.0). Each point 
has been obtained by fitting the numerically projected Einasto profile with 
UA\\ for b = 2n and Ci = 1. 



Cl = l 
6 = 2n 

X = R/r^2 

and, C,2 and /i are numerically estimated (Fig.l and Fig.2) for the 
specific conditions of 6 = 2n and (^i = 1 to yield: 

0.05657 0.00903 



C2 



1.1513 



1.5096 0.82505 
+ :r- 



0.66299 



(2.12) 



(2.13) 



We would also like to note that based on the accuracy of an approx- 
imation needed, one can neglect the contribution from region III 
(the 5 term) and with ("i = 1 fit only for ^2- Although we have not 
explored a functional relation of C,2{n), we note that the best-fitting 
^2 varies weakly from 1.2145 at n = 10 to 1.2194 at n = 3.0 to 
1.2424 at n = 1.0, with a maximum error within 30 r_2 (usually 
reaching a peak around r2oo) of 0.6% at n = 10, 1.2% at = 3 
and 2.2% at n = 1. Hence, within the current domain of 3 < n < 8 
from N-body simulations (corresponding to 0.12 < a < 0.3) one 
can neglect the 5 term and set (^i = 1, and C,2 ~ 1.2176 at the cost 
of an error in the range (0.5% to 1.5%). 

The parametrizations ( |2.12t and J2.13l > worked very well even 
at n=14 and n=0.95. However, we have not tested n > 10 and n < 
1 values rigorously. If greater accuracy is needed, we recommend 
fitting for (2 and /.i. [Refer to the discussion following l |2.9l ), on why 
the S term is negative]. This is a useful result because the surface 
mass density is expressed entirely as a function of the 3D spatial 
density parameters, which is not the case for any existing projected 
fitting function. Equation ( |2.10t or J2.1 U thus serves two purposes: 
One, given a 3D Einasto profile, J2.10t gives a good approximation 
to its 2D surface mass density. And two, a good fit to some 2D 
observations with ( 12. lot for example surface mass density from 
lensing, will give us the 3D spatial density parameters of a Einasto- 
like profile. 

This is also quite unlike fitting 2D Einasto profiles with Sersic- 
like functions, where one first needs to fit for the 2-D and then de- 
project to fit for the 3D shape parameters (M05), which usually are 
different without any known existing functional relation between 
them. 



3 COMPARISON WITH NUMERICALLY PROJECTED 
EINASTO PROFILES 

3.1 Estimating (2 and /i for Eb 

We numerically integrate l l2.lt for the profile in il.H and obtain 
Ejv(i?) in the domain R : {0 — 30) r_2 for 90 profiles with a 
shape parameter in the range 0.1 ^ a ^ 1. A resolution in R of 
0.002 r_2 (~ 0.05 rconv), allows us to quantify errors due to our 
approximation in a domain R « rconv, and we report compari- 
son of errors up to 30 r_2 or up to a R where Ejv(-R) ~ 10~*E(0) 
whichever is earlier. 

With b — 2n and (^1 = 1, we use a non-linear least squares 
Levenberg-Marquardt algorithm to estimate the best-fitting values 
for ^2 and /i by fitting i2.11t to each of the numerically generated 
Ejv profiles. We find that (2 l l2.12t and fj, l |2.13t are best described 
by a second and third degree polynomial respectively, in a = 1/n 
(Fig. 1 and Fig.2). 



3.2 A detailed look at errors in the E e approximation 

Fig. 3 and Fig. 4 describe the fractional error profile between our 
model Ee and the numerically projected Einasto profile Ejv(R) 
for a wide range of the shape parameter a = 1/n. Fig. 3 is more 
relevant to current N-body simulations, where in a seems to be in 
the range 0.1 to 0.25. It is worth noting here, that for a as high as 
0.25 (n = 4.00), the largest errors are < 0.3% in the range (0 to 
30) r_2. 

In case the range 1 n ^5 3 becomes relevant in the future, 
where N-body haloes for a > 0.4 have not yet been found, we also 
present in Fig. 4 a comparison between E_e and Ejv. This is also 
the domain where our assumptions in the z > R region are weaker. 
Nevertheless, the accuracy of the approximation is striking with the 
worst error < 1% within R < 5 r_2. 



3.3 Comparison with Sersic profiles 

In this section we discuss results of fits to Ejv with a Sersic func- 
tion and superimpose the E e model (black solid lines) I I2.1U for a 
relative comparison between Eg and E_b (Fig. 5 and Fig. 6). The 
fits presented here were obtained using log scale of density. One 
can also obtain fits with density in linear scale; a comparison of 
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Figure 3. Fractional Error between and S]v for 0.1 ^ o 0.25. 
4 ^ n ^ 10. The x-axis plotted in log scale is expressed as a ratio of the 
2D projected radius R and the 3D scale radius r_2 of the corresponding 
Einasto profile. r200 ~ 6 r_2 while rconv ~ 0.05 r_2. 
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Figure 4. Fractional Error between S^; and Sjv for 0.33 ^ a ^ 1.00, 
1 ^ n ^ 3. The x-axis plotted in log scale is expressed as a ratio of the 
2D projected radius R and the 3D scale radius r_2 of the corresponding 
Einasto profile. r200 ^ 6 r_2 while rconv ~ 0.05 r_2. 



residuals for the case a — 0.17 is presented in Fig. 7. The best 
fit parameters and consequently the error profile are quite differ- 
ent and is an indication that the best-fitting Sersic profile does not 
provide an adequate representation of a projected Einasto profile. 

Not evident in the plots (Fig. 5-7) are the extremely large errors 
in the central density Eg (0) when fitting using log density. For a = 
0.1 (Fig. 5) the relative error at 7? = is 375%, for a = 0. 17 (Fig.5 
and 7), the relative error is 215% and 33% for a = 1. However, we 
observed empirically that fits in log density happen to reflect the 
significance of Re as enclosing half the total mass of Sat (it is to 
be noted that Re by definition encloses half the total mass of Eg). 
This is because the domain over which Eg (with fits obtained in log 
density) overestimates the density is a relatively small contributor 
to the total mass and fits with log density in the region R > r_2 
are good. Since, we can reliably estimate one of the Sersic profile 
parameters (Re) through fits with log density (as opposed to none 
in linear scale), we have presented results of fits in log density in 
Fig. 5 and Fig.6. 

The large errors in the central density also have serious con- 



sequences for strong lensing. Typically the strong lensing regime 
extends up to ~ 0.1 r_2. In strong lensing, image positions cor- 
respond to extremum (minima, maxima and saddle) points of the 



time delay surface. The j 
by: 

dij ~ Pi + ■ 



image 9ij for the i source Pi is given 



i9ij 



(3.1) 



where k{9) is the normalized surface mass density at an angular 
position 0. 

Consequently contributions to the integral in l l3.lt from un- 
usually large density near the center (as a result of fits with a Sersic 
profile in log density), will produce image separations larger than 
what one can expect from a numerical projection of the 3D Einasto 
profile. 

Sersic profile fits to model a projected Einasto profile should 
thus be avoided especially if the central region (i?<0.1r_2) is be- 
ing excluded. This is because, as shown in Fig. 7, one can get a 
reasonably good fit (relative error within 10%) over a large range 
of 7?>0.1r_2 (fitting in log density) giving an indication that the 
Sersic profile is a good representation of projected Einasto profile, 
but doing so will lead to even larger errors in the excluded central 
region (i?<0.1r_2). One should thus use caution in interpreting 
the other structural parameters, the shape parameter m or A and the 
central density Eg(0). 

Fits in linear density present a different problem. Even though 
the central errors are much better (relative error ~10% for a=0.17. 
Fig. 7) than fits using log density, the best-fitting Re does not en- 
close half the total mass of Ejv and the errors for large R keep 
increasing with R. 

Although the Sersic profile is not a good representation of a 
projected Einasto profile, one can fit Sersic profiles in limited do- 
mains of R and obtain an estimate of the shape of the projected 
Einasto profiles in those domains only, but be careful to not use the 
resulting best-fitting values of Re and Eg(0) as a true represen- 
tation of the half-mass radius and central density of the projected 
Einasto. In Fig.8, we present two such relations between the 2D 
Sersic index (A= — ) and the 3D Einasto index a=- in two do- 
mains R< r-2 and R> r_2. In these 2 regions, A(a) can be de- 
scribed as power laws. For the domain R< r_2 we find: 



A(a) 



1.332 a 



and for R > r_2 defined as in §3.1: 

A(q) = 1.037 a° 



(3.2) 



(3.3) 



Not only is the shape (A) different in the two domains, their evolu- 
tion with a is also different. Nevertheless, this result can be useful 
in obtaining an estimate of the shape of projected Einasto profile in 
these two domains demarcated by r-2. 

The E_E model, with errors < 0.5% does not face any of the 
above issues. Further, unlike the Sersic profile, the Eb model can 
predict the central density with almost 0% errors due to the exis- 
tence of an analytical solution. Thus, if the underlying 3D distri- 
bution is Einasto-like, the 2D distribution should be modeled with 
Eb. 



3.4 Extracting the 3D-paraineters from the E e model 

The Eb model is not just a good description of the projected 
Einasto profile but is also expressed in terms of the 3D Einasto 
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Figure 5. Fractional Error between the best-fitting Sg and Sjy witii 
Einasto index 0.1 ^ a ^ 0.25, 4 ^5 n ^ 10 with fractional error for 
the best-fitting S^; (black solid Hnes 0) superimposed. The x-axis plot- 
ted in log scale is expressed as a ratio of the 2D projected radius R and the 
3D scale radius r_2 of the corresponding Einasto profile. r2oo ~ 6 r_2 
while Tconv ^ 0.05 r_2. 




Figure 6. Fractional Error between the best-fitting Eg and Sjv with 
Einasto index 0.33 ^ a ^ 1.00, 1 ^ n ^ 3 with fractional error for 
best-fitting S^; (black solid lines fa 0) superimposed. The x-axis plotted in 
log scale is expressed as a ratio of the 2D projected radius R and the 3D 
scale radius ■r_2 of the con'esponding Einasto profile. r2oo ~ 6 r_2 while 

Tconv ~ 0.05 r_2. 



profile parameters. It should thus be possible to recover the 3D pa- 
rameters (q, r_2, p-2) from fits to 2D distributions that subscribe 
to an underlying 3D Einasto-like system. 

For the wide family of numerically projected Einasto profiles 
Ejv described in this paper, we could recover the 3D parameters 
for all of them with an accuracy of ~ 10"'' or better, by fitting E n 
with J2.1U and the parametrizations of ( I2.12t and J2.13b through 
a non-linear least squares Levenberg-Marquardt algorithm. Given 
that, as of now robust data (within virialized regions) from N-body 
simulations are in the domain Vconv to r2oo, the fits were performed 
in this domain. In passing, we note that our results are even better 
if we fit from to 30 r_2- Such a high degree of accuracy indicates 
that if the 2D distribution is indeed a projected Einasto profile, the 
3D parameters can be recovered very well even from a limited ra- 
dial range of observations. 

We note that, for the entire range of 0.1 ^ {a — i) ^ 1 the 



Figure 7. Fractional Error in S5 from fits to Ejv (for = 0.17) with a 
Sersic profile using log density (dot-dashed line) and a Sersic profile using 
linear density (dashed line). Not shown in the plots are the errors in cen- 
— . The errors are 215% for the fit with log density and 



tral density 



10% for the fit with linear density. The fractional error in Eg (black solid 
line 0) is superimposed for comparison. The x-axis, plotted in log scale, 
is a ratio of the 2D projected radius R to the 3D scale radius r_2 of the 
corresponding Einasto profile. r2oo ~ 6 r_2 while Tconv ~ 0.05 r_2. 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Einasto (3D) shape parameter a=1/n 

Figure 8. Sersic index(A=l/m) of a projected (2D) Einasto profile as a func- 
tion of the 3D Einasto index (a=l/n) in two domains R < r—2 (dotted) and 
ij > r_2 (solid) indicating the absence of a unique Sersic index A for all R. 
The power law relations for \(a) in equations j3.2) and j3.3) are applicable 
only in these domains. 



fits always converged for an intial guess of a < true a, rs > true 
r_2 and p_2 in the range (0.1 to 4) An inital guess of very 

low a ~ 0.05 and a guess for r_2 ~ r2oo for the type of object 
(galaxy or cluster) being considered can be a reasonable starting 
value for the fit to converge. We also did not encounter any local 
minima, i.e. if the fit converges, it always converged to the true set 
of (a, r_2, P-2)- 



4 SUMMARY AND CONCLUSIONS 

Non-parameteric estimates of density profiles in N-body simula- 
tions (Nav04,M06) favour Einasto-like profiles, since they pro- 
vide better fits than the two-parameter NEW and Moore profiles. 
iMerritt et alj ( l200^ have also shown that a de-projected Sersic pro- 
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file fits the 3D halo mass distribution almost as well as the Einasto 
profile, and a Sersic profile provides good fits to non-parametric 
estimates of surface mass densities (M05) of the Nav04 N-body 
haloes. 

We have observed that fits with a Sersic function (Eg) to a 
numerically projected Einasto profile (Eat) are sensitive to whether 
one fits using linear density Es (errors increasing for large R) or 
log density /n(Es) (errors increasing for small R) yielding widely 
varying results. Consequently, the Sersic profile does not give an 
adequate description of the projected Einasto profile. 

Sersic profile fits to the surface mass density of N-body haloes 
(M05), whose 3D spatial densities are well fit by Einasto profiles 
with 0.12^0^0.22, have been obtained from the limited radial 
range of rconv to r2oo. For the haloes in M05 and Nav04, this range 
is generally less than two decades in radius. Hence, if the 3D dis- 
tribution is indeed Einasto-like, interpreting structural properties 
from fits with a Sersic profile, especially m in ( ll.2t as the shape 
parameter, Es(0) as the central density and Re as the half-mass 
radius, can be misleading. 

In this paper, we have provided an analytical approximation 
1 I2.I It to the projected surface mass density of Einasto-like 3D 
density distributions. The fit errors are well contained to <2% for 
the projected radial range O^Ci?^ (10 — 30) r_2 equivalent to 
(3 — 5)r2oo and shape parameter, l^Cn^ClO, or O.l^a^Cl. This 
model can therefore be used both as a fitting function for 2D ob- 
servations and also to extract the 3D parameters of Einasto-like 
profiles. Since E_b fits a projected Einasto profile in a wide radial 
range, it can be used for fitting strong and weak lensing observa- 
tions in systems whose total 3D density distribution is believed to 
be Einasto-like. One can also numerically integrate l l2.11b to get 
reliable estimates of the mass enclosed. 

Finally, we note that the form similarity of dl.ll l and l |1.2| >. i.e. 
fitting functions that describe the 3D mass density of dark matter 
haloes and the 2D light distributions of galaxies, respectively, could 
be largely coincidental and should be used with caution when draw- 
ing conclusions about the similarity of dynamical evolution that 
lead to the formation of the stellar components of ellipticals and 
dark matter haloes. 
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